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Abstract 

Radial velocity dispersion in spheroidal and elliptical galaxies, as a 
function of radial distance from the center of the galaxy, has been derived 
from Cosmological Special Relativity. For velocity dispersions in the outer 
regions of spherical galaxies, the dynamical mass calculated for a galaxy 
using Carmelian theory may be 10 to 100 times less than that calculated 
from standard Newtonian physics. This means there is no need to include 
halo dark matter. The velocity dispersion is found to be approximately 
constant across the galaxy after falling from an initial high value at the 
center. 

Key words: Cosmological General Relativity, elliptical and spheroidal galax- 
ies, dispersion velocity 

1 Introduction 

The motion of stars or the motion of gases as characterized by the spectro- 
scopic detection of neutral hydrogen and other gases in spheroidal and elliptical 
galaxies has caused concern for astronomers for many decades. Newton's law 
of gravitation predicts much lower radial velocity dispersions than those mea- 
sured. This has led to the assumption of the existence of halo 'dark matter' 
surrounding galaxies but transparent to all forms of electromagnetic radiation. 

Carmeli [31 Sj formulated an extension of Einstein's general theory, in an 
expanding universe taking into account the Hubble expansion as a fundamental 
axiom, which imposes an additional constraint on the dynamics of particles [T]. 
Carmeli believes the usual assumptions in deriving Newton's gravitational force 
law from general relativity are insufficient, that gases and stars in galaxies are 
not immune from Hubble flow. As a consequence a universal constant ao is 
introduced as a characteristic acceleration in the cosmos. 

Using spherical coordinates Carmeli [2] successfully provided a theoretical 
description of the Tully-Fisher law for spiral galaxies. Following Carmeli's lead, 
Hartnett [5] showed that the same line of reasoning leads to plausible galaxy 
rotation curves in spiral galaxies using cylindrical coordinates. 
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In this paper I take the analysis further, using spherical coordinates and 
model the gravitational potential and the resulting forces determining how test 
particles move in spheroidal and elliptical galaxies with an appropriate density 
distribution that reflects the observed luminous matter distribution. Similar to 
what was found in spiral galaxies , two acceleration regimes are apparent here 
also. In one, normal Newtonian gravitation applies. In that regime the effect of 
the Hubble expansion is not observed or is extremely weak. It is as if particle 
accelerations are so great that they slip across the expanding space. In the other, 
new physics is needed. There the Carmelian metric provides it. In this regime 
the accelerations of particles are so weak that their motions are dominated by 
the Hubble expansion and as a result particles move under the combined effect 
of both the Newtonian force and a post-Newtonian contribution. 



2 Gravitational potential 

In the weak gravitational limit, where Newtonian gravitation applies, it is suf- 
ficient to assume the Carmeli metric with non-zero elements goo = 1 + 
(744 = 1 + gkk = —1, (fc = 172,3) in the lowest approximations in both 

1/c and 1/r. Here a new constant, called the Hubble-Carmeli constant, is in- 
troduced T « I /Ho- See Ref. [Ij for details. The potential functions $ and ^ 
are determined by Einstein's field equations and from their respective Poisson 
equations, 

= 47rG/9, (1) 
V^* = (2) 



where p is the mass density, oq a universal characteristic acceleration Oq = c/r 
and c is the speed of light in vacuo. 

A comparison of $ and ^E* in ([T]) and ^ leads to ^E* — ^/flg within an 
arbitrary additive constant. Since both potentials are defined with respect to 
the same co-ordinate system, in reality, we only need deal with one potential 
function, the gravitational potential, <&. 

The density function that best describes the luminous distribution of matter 
in spherical and elliptical galaxies is 

M r, 

= ~A 2l~^ ^2' 

47r r'^(r + rc) 

where Vc is the core radius which contains half of the luminous matter. The 
parameter M is the total mass of the galaxy as measured at r = cx). Therefore 
the mass at radius r is 

Spherical symmetry has been assumed to simplify the problem. To consider 
fully oblate spheroids would introduce some asymmetry to the problem but in 
principle the solution found here should broadly apply to any spherical matter 
distribution and would also apply, to first order, to oblate spheroids. 
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In spherical coordinates (r, 0, the potential $ that satisfies ([T]) can be 
found by integrating over ^ using ([3]). Only the radial dependence remains. 

, GA/ r, /l + re/A\ re/A 1 



1 + rjr J 1 + rJA 

where A is the radial extent of the matter distribution. This may be approxi- 
mated where A ^ rg, which is the usual case. Hence we get, 

^r) = -—log(l+^-^). (6) 



3 Equations of Motion 

The 5D line element for any two points in the CGR theory is ds^ = gooc^dt^ + 
gkk{dx^Y' + gi/^T^ dv^ , where k = 1,2,3. The relative separation in 3 spatial 
coordinates — {x^Y + {x^Y + {x^Y and the relative velocity between points 
connected by ds is v. The Hubble-Carmeli constant, r, is a universal constant 
for all observers. The equations (B.62a) and (B.63a) from Ref. [1] are used in the 
sequel to derive the appropriate equations of motion to lowest approximation 
in 1/c. 



3.1 Newtonian 

It follows from Carmeli's equation (B62a) using ([3]) and (O, and the usual form 
of the circular motion equation 



v"^ d<I> 
r dr ' 

that 



(7) 



2 GM 1 

r 1 + Tc/r ^ ' 

where G denotes the gravitational constant. 

Equation ([S]) is the usual Newtonian result for the speed of circular motion 
in a spherical gravitational potential. This equation has been plotted in curve 
1 of fig. mja) as a function of radial position from the center of a galaxy in kilo- 
parsecs {kpc) where kpc « 3.08 x 10^^ m. Figure [ijb) shows the corresponding 
accelerations. Curve 1 is the Newtonian acceleration. Through out this paper 
M is expressed in solar mass units Mq « 2 x 10"^" kg. 



3.2 Carmelian 

Using '5 = ^/oq in Carmeli's (B63a) results in a new equation 

r dr 
Jo 

which has been integrated and solved for w as a function of r. Using the potential 
$, determined from ^ and in results in 

2 r3/2 
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Figure 1: (a) Orbital velocity of a test particle in a dwarf spheroidal galaxy 
calculated from jS]) and ^ and where M = IO^A/q and = 0.004 kpc, 
indicated by the black vertical line; (b) Corresponding acceleration of a test 
particle for both the Newtonian and Carmelian regimes. Note they coincide at 
the critical acceleration 2/3ao. 



4 



which describes the expansion of space within a galaxy. 

In Refs [2] and [5], using spherical co-ordinates, it was found that in the limit 
of large r and where all the matter was interior to the position of a test particle, 
such a particle is also subject to an additional circular motion described by (fTU|) . 
Apparently this is the result of the expansion of space itself within the galaxy 
but in an azimuthal direction to the usual center of coordinates of the galaxy. 
In this paper also the same result PH)) was obtained. 

Carmeli [3] determined a TuUy-Fisher type relation using the Newtonian 
circular velocity equation expressed in spherical coordinates, 

^ ^. (11) 

r 

where it is assumed that test particles orbit at radius r outside of a fixed mass 
M. Then by eliminating r between ([Tl]) and (fTO|) we get the result. This is 
achieved by taking the 3/2 power of pTjl and multiplying it with (|10p yielding 

= GM^ao. (12) 

So by applying the same approach with ([5]), raising it to the 3/2 power 
and multiplying it with (fTU|) . we can derive an equation describing the circular 
motion of test particles in spherical galaxies. The result is 

But in elliptical and spherical galaxies we generally don't observe any group 
circular motion. Therefore the trajectories of the individual orbits, though 
circular around the center of mass, are randomly oriented. Nevertheless we can 
calculate the orbital speed of a test particle from Curve 2 of fig. [H^a) 

shows one for the same parameters as is assumed for the Newtonian curve, as a 
function of radial position from the center of a galaxy. However for an individual 
particle if r 3> Tc equation ()13p recovers the form of the TuUy-Fisher relation 

m- 

The corresponding acceleration may also be calculated from w^/r using (fTH]) 
and is shown in curve 2 of fig. [Ijb) for the same galaxy parameters. It is also 
compared with the critical acceleration 2/3ao. 

4 Radial velocity dispersion 

Assuming an spherically isotropic dynamic cloud of gases and stars, where the 
radial velocity dispersion is equal to the or ip dispersion (i.e. — ^) the 
hydrostatic Jeans equation becomes 

±{pir)a',) = -pir)^. (14) 
ar ar 

4.1 Newtonian 

Using ([7|), ([H]) and ^ we solve the differential equation (fTi|) for the radial 
velocity dispersion 

, lGAf„, , , , 

'^r = 2(x), (15) 

Z Tr 
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Figure 2: Dimensionless velocity dispersion squared. Curve 1 is S(x) from (fTC|) 
and curve 2 is r(a::) from ([20|) as functions of the relative radius x = r/vc- 



where 



E{x) = (1 + 2a;)(l - 6a;(l + x)) + 12x^(1 + xflog 
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(16) 



and X = t/tc- The dimensionless function S(r/rc) is shown as curve 1 in fig. [21 
For r 3> the function S(a;) w 2/52; which means (jl5p becomes 

, GM 



5r 



(17) 



This is the same result determined from the virial theorem. 

So in the limit where r —^ oo, the function 2 — > and where r —^ 0, the 
function S ^ 1. From Newtonian theory we expect that the radial velocity 
dispersion in spheroidal and elliptical galaxies to tend to zero where the gravi- 
tational acceleration is very weak. And where the gravitational acceleration is 
strongest (r = 0) the radial dispersion is maximum and equal to the rotation 
velocity of a particle at twice the core radius. 



4.2 Carmelian 

Using (jl3p we can construct an equivalent gravitational potential which includes 
the effects of the Hubble flow on test particles in the galaxy. By taking u^/r 
and using this in (fTi]) instead of the gradient of the Newtonian potential we get 

1/2 



9 1 / 2 



T{x), 



(18) 
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where 



T{x) 
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:3/^(l + a;)i/*{7 - 4x(21 - 32(a;(l + xfY^^ + ... 

... + 8x(7 + 4T-4(x(l + a;)3)i/*))}, 



(19) 



and X — r/vc- See curve 2 in fig. [2] for T{r/rc). 

In the limit where r ^ oo, the function F ^ 1 and where r ^ 0, the function 
r — !■ 0. From Carmehan theory then we expect that the radial dispersion in 
spherical and elliptical galaxies to tend to zero where the gravitational accel- 
eration is strongest (r = 0) and to tend to a constant where the gravitational 
acceleration is very weak. It also follows from (fT8| that for r ^ Tc 



4.3 Discussion 

Similarly to what was seen in the case of rotation curves fS], two regimes must 
be considered. When accelerations are greater than the critical acceleration 
2/3ao we have the Newtonian regime and we expect the velocity dispersion to 
follow (|15p but when accelerations are less than 2/3ao we have the Carmelian 
regime and we expect the velocity dispersion to follow (fT5)) . Notice that even 
though the velocity dispersion from the Newtonian calculation p^]) scales with 
core radius (re), the velocity dispersion from the Carmelian calculation (llSp is 
independent of core radius. 

The Newtonian regime applies where the Newtonian acceleration qn > 
2/3ao and the Carmelian regime applies where the Carmelian acceleration ac < 
2/3ao. In fig. [3]radial velocity dispersion curves are shown for a dwarf spheroidal 
(curve 1) with M = 10^ Mq and a massive elliptical (curve 2) with M — 
10^^ Mq. Only the regions of the curves which represent valid velocity dis- 
persion are retained. The results indicate that, as a function of radius from 
the center of the galaxy, the velocity dispersion should fall from a central high 
value to a constant value as r becomes much greater than rc. This is commonly 
observed in these type of galaxies. 

4.3.1 Dwarf spheroidal galaxies 

It follows from (PO)) . that we can write an equation for the radial velocity dis- 
persion that we expect to observe in a typical dwarf galaxy, of mass M (in units 
of Mq ) as 



assuming the observations are obtained from the sources far from the core radius. 
That is, sources that are in the Carmelian acceleration regime. Here a value of 
T = 4.276 X 10^^ s has been assumed p^, which is equivalent to — 72.17 
km s^^ Mpc^^ , and hence qq — 4.674 x 10"^" m/s . 

By calculating the mass required for typically observed values of cr,. , in dwarf 
spheroidal galaxies, according to the both Newtonian (fT5|) and Carmelian equa- 
tions (fT8|) developed here a comparison can be made. The galaxy mass (M) 
calculated for different radial velocity dispersions <Tr at x — r/vc — 10 (for a 
typical core radius of Tc = 0.004 kpc) are listed in Table I. 
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Figure 3: Left axis: Radial velocity dispersion as a function of x = r/rc- Right 
axis: Temperature of gas in units of keV. Curve 1 is ar{x) for a dwarf spheroidal 
galaxy with M = 10^ Mq and curve 2 is ar{x) for a massive elliptical galaxy 
with M = 10^3 



Table I: Masses of spheroidal dwarf galaxies from 
velocity dispersion cr^ at a; = r/rc = 10 



(7r 


Carmelian 


Newtonian 




[km/s] 


M[x10^Mq] 


M[x10^Mq] 


Ratio 


4.40 


1.00 


9.38 


9.38 


6.58 


5.00 


20.9 


4.20 


7.82 


10.0 


29.7 


2.96 


11.7 


50.0 


66.2 


1.33 
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4.3.2 Bright central elliptical galaxies 

Similarly from (|20p we can write an equation for for the radial velocity dispersion 
that we expect to observe in a typical massive elliptical galaxy, of mass M (in 
units of Mq) as 

= 247.4 1^^^ j (22) 

also assuming the observations are obtained from the sources that are far from 
the core radius, in the Carmelian acceleration regime. 

Again by calculating the mass required for typically observed values of cr^, 
in bright ellipticals using both models a comparison is made results listed in 
Table II. In this case a typical core radius of = 40 kpc was assumed. 



Table II: Masses of bright elliptical galaxies from 
velocity dispersion ar Sii x — r/vc — 10 





Carmelian 


Newtonian 




[km/s] 






Ratio 


139.1 


0.1 


9.38 


93.8 


247.4 


1.0 


29.7 


29.7 


369.9 


5.0 


66.3 


13.3 


439.9 


10.0 


93.8 


9.38 



4.3.3 Temperature 

By considering the random motion of gas particles we can write an expression 
for the temperature of a cloud of gas in the galaxies that is heated by the 
gravitational potential the gas feels whether it be in the Newtonian or Carmelian 
regime. Using /i = 0.609 as the mean atomic weight of the gas we can equate 

/impCT^ = kT, (23) 

where k is Boltzmann's constant, Wp is the proton mass and T is the temper- 
ature of the gas. Using (|23p the temperature of the gas in a galaxy heated 
in this manner can be determined from the radial velocity dispersion. This is 
shown on the right axis of fig. [H The result indicates that the temperature 
of the gas across spherical galaxies is expected to be approximately isother- 
mal. Also, due to the thermal heating of intra-galactic gas, massive elliptical 
and spheroidal galaxies emit radiation in the X-ray part of the spectrum, while 
dwarf spheroidals emit in the near infrared. This is consistent with observations. 

5 Conclusion 

This theory suggests that it is the Carmelian regime that is applicable at low 
accelerations where r 3> and hence the masses of galaxies are overestimated 
from the observed dynamics. In Table I and II typical dispersion velocities are 
used but the core radius is fixed at = 0.004 kpc and Tc = 40 kpc, respectively. 
So depending on the exact details for a galaxy one may get 10 to 100 times 
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more mass from a Newtonian calculation than from a Carmelian calculation. 

This brings the masses more in hne with estimates from the himinous material. 
Therefore the need to invoke halo dark matter is avoided. It is also true that ct^ 
is observed to fall from a central higher value near the center of the spheroidal 
galaxy and become constant as a function of radius from the center. This is 
indicated by a Newtonian regime in the center which becomes Carmelian as 
a function of radius. In massive ellipticals radial velocity dispersion is often 
observed to be approximately constant as a function of radius. This is indicated 
by observations where most of the galaxy is dominated by the Carmelian regime. 
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